Coarsening and Slow-Dynamics in Granular Compaction 
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We address the problem of the microscopic reorganization of a granular medium under a compaction 
process in the framework of Tetris-like models. We point out the existence of regions of spatial 
organization which we call domains, and study their time evolution. It turns out that after an initial 
transient, most of the activity of the system is concentrated on the boundaries between domains. 
One can then describe the compaction phenomenon as a coarsening process for the domains, and 
a progressive reduction of domain boundaries. We discuss the link between the coarsening process 
and the slow dynamics in the framework of a model of active walkers on active substrates. 
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The phenomenon of granular compaction involves the 
increase of the density of a granular medium ^ sub- 
ject to shaking or tapping. Triggered by experimental 
results of the Chicago group that suggested that 
compaction follows an inverse logarithmic law with the 
tapping number, several models have been proposed to 
explain the slow relaxation features of granular media 
|^-|lO|. Though in all these different cases a very slow 
relaxation (eventually logarithmic) is reproduced, an ex- 
plicit connection between the above models and a real 
granular medium is however still rather tenuous. 

The aim of this paper is to elucidate the origin of the 
very slow relaxation studying explicitly the microscopic 
response of a granular medium subject to shaking. We 
address this problem within the framework of the recently 
introduced class of Tetris-like models |^ which are known 
to reproduce several features observed experimentally in 
granular materials such as slow dynamics, segregation, 
aging and hysteresis. 

We find, quite surprisingly, that the system reorganizes 
under the shaking dynamics into several ordered regions 
(see [|ll| , |l2| for other examples). We call these domains, 
and study their time evolution. After a short transient, 
most of the activity of the system is concentrated along 
the boundaries between domains (we note that this con- 
cerns only a small part of the entire system). Under 
shaking, the domain boundaries move throughout the 
system and free the vacancies they encounter leading to 
a progressive densification. Moreover, when two domain 
boundaries meet, they annihilate. One can thus describe 
the compaction phenomenon in this system, as a coars- 
ening jl^ (i.e. a domain growth) process. As the system 
compactifics, the domains coarsen and the boundary re- 
gions are reduced and thus the process becomes slower. 
We give a quantitative description of this phenomenon 
studying the behavior of the space-time correlation func- 
tion that is expressed by C{r,t) ~ f{f/£,{t)) where ^(t) 



is the correlation length, say the typical size of a do- 
main. This coarsening of domains is related to the slow 
compaction process by measuring the persistence [|4j ex- 
ponent of the phenomenon, as well as by measuring the 
activity and the motion of the domain boundaries. 

Let us briefly recall the definition of the Tetris model. 
Although this class of models allows for an infinity of par- 
ticles types, shapes and sizes, here we use, without loss 
of generality for the main features, a system of elongated 
particles. These occupy the sites of a square lattice tilted 
by 45°, with periodic boundary conditions in the horizon- 
tal direction (cylindrical geometry) and a rigid plane at 
its bottom. Particles cannot overlap and this condition 
produces strong constraints (frustration) on their relative 
positions. The system is initialized by inserting grains 
at the top of the system, one at the time, and letting 
them fall down, performing under the effect of gravity, 
an oriented random walk on the lattice, until they reach 
a stable position, i.e. a position from which they cannot 
fall further. The effect of vibrations is implemented by 
means of a two-steps Monte-Carlo algorithm mimicking 
a tapping procedure. The role of the tapping amplitude 
is played by a parameter < x < 1 that describes the 
strength of the bias in the particle movement, induced 
by the gravity (we refer to for the details). 
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FIG. 1. Time-evolution of the Tetris model with x = 0.5 
resolved in antiferromagnetic domains at times (number of 
taps) t — Q, 10^ and lO'^ from top to bottom. White squares 
are voids while red and blue squares represent particles be- 
longing to the two possible domain types. 
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In the simplest version, the Tetris model consists of 
a single rod-like type of particle (rectangles of uniform 
size a X b with a ~ 0.75 mesh units and 6 = 0) with 
two possible orientations (along the principal axis of the 
lattice) chosen to be equally probable. A generic config- 
uration can be described by assigning to each site of the 
lattice (x, y) {0 < x < is the horizontal coordinate 
and < y < Ly is the vertical one), a variable a{x,y,t) 
whose value is if the site is empty and ±1 if the site is 
occupied by a particle with one of the two possible orien- 
tations. At every density the system can be resolved into 
domains, i.e. regions in which the staggered magnetiza- 
tion keeps a definite sign and each domain presents an 
antiferromagnetic order with vacancies (-1-1 particles on 
odd (even) rows and —1 particles on even (odd) rows). 
One can then observe the evolution of the compaction 
dynamics in terms of the evolution of the domains. 

At the beginning, after pouring the grains into the con- 
tainer, i.e. in the so-called loose density state, the system 
presents a disordered structure with an alternation of the 
two types of domains, even though the domain bound- 
aries do not as yet span the system from top to bottom. 
At this stage the number of domain boundaries depends 
on the aspect ratio (height by width) of the container: the 
smaller the aspect ratio (wider is the system) the larger 
the number of domains. The domain size is of the order 
of the height of the system, and is almost independent of 
the horizontal size of the system L^, as long > Ly. 

The compaction can now be seen as a slow elimina- 
tion of the voids frozen in the different domains. Since 
the system changes only at the domain boundaries, a 
void in the bulk of the domain can be freed only when it 
comes in contact with a domain boundary. The domain 
boundaries are then the only regions where the activity 
of the system is concentrated. Fig.|l| shows an example of 
the time-evolution of the Tetris model resolved in anti- 
ferromagnetic domains jisll . It is important to note how 
narrow systems {L^ <^ Ly) may display a pathological 
behavior (blocking) if the system has an almost single- 
domain like packing. 

Let us now describe the coarsening dynamics in a more 
quantitative way. We have monitored the evolution of the 
(longitudinal) correlation function defined as: 

P y = X = 

where Np is the number of particles in the bottom half of 
the system. A pair of particles inside the same domain 
gives a positive contribution to C(r, t) while a pair in 
different domains gives a negative contribution. With 
this definition, the correlation function is not sensitive 
to density changes and reflects only the evolution of the 
domain sizes. 

We perform extensive simulations of wide systems (to 
avoid blocking) and we try a scaling collapse of several 



C{r,t) curves at different times: i.e. we look for a char- 
acteristic length ^{t) such that C{r,t) = f{r/^{t)), the 
length ^ (t) representing the average (horizontal) domain 
size. As in standard coarsening dynamics, ^{t) grows in 
time and when ^(t) w L^ the growth stops (blocking for 
single domain systems). 
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FIG. 2. Scaling collapse of (1) for a system with Lx = 800 
and Ly — 50, simulated up to a time t — 10^ with x — 0.1, 
and averaged over 200 different dynamics, ^{t) is chosen such 
that C{^{t),t) — Co for all the curves obtained at different 
times for Co = 0.2 (One gets the same result for a large range 
of values of Co). In the inset is shown ^(t) vs. t for differ- 
ent values of a; = 0.1,0.5,0.9. The crossover occurs at an 
a;-dependent t* such that (,{t*) ~ Ly. In the second regime 
one observes an exponent equal to 0.5 only for a; = 0.9 while 
slightly smaller exponents are observed for smaller x. We 
believe these deviations from an exponent 0.5 are transients 
evolving towards the asymptotic diffusive value. 

However, the height, Ly, of the system is another char- 
acteristic scale, and depending on whether ^ is smaller 
or larger than Ly, two different regimes can be observed. 
A quantitative analysis of ^ versus time reveals that for 
C(t) < Ly, ^{t) oc whereas for ({t) > Ly, we observe 
a faster growth : £,{t) oc t^-^ (See Fig.||). 

These results can be interpreted according to the fol- 
lowing scenario. One can imagine in general that the 
size S{t) (area) of the domains grows as a power of time. 
The correlation length ^(t) is a measure of the lateral 
size of the domains. Now, as long as the domains grow 
in an isotropic way we can expect the area to scale as 
S{t) ~ £,{t)'^- In this regime the domain walls, though 
biased by gravity, do not yet span the system in the y 
direction. Later on, when the domain walls span the 
system from top to bottom, the coarsening dynamics is 
dominated by the diffusion of these almost vertical walls, 
which eventually collide and annihilate each other. At 
this stage we expect the area scales as S{t) ^ Ly£^{t). The 
crossover we observe is then compatible with a growth 
S{t) ~ which gives $,{t) ~ t^^"^ in the early regime 
followed by the asymptotic behaviour ^(i) ~ t^^"^. The 
crossover between the first and the second regime is evi- 
dent in the first two pictures of Fig.0. 
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FIG. 3. R{t) (see text for the definition) vs. time for difler- 
ent values of x. The inset sliows tlie local logaritiimic slopes. 

How now do we relate this coarsening behaviour to 
the density relaxation? The motion of the domain walls 
occurs in a background of vacancies. We observe that 
regions swept through by the domain walls compactify 
by triggering particle rearrangements, while regions not 
yet swept through remain disordered. Thus within each 
domain there is a compactified region through which the 
domain wall has made several forays, and a disordered 
loose region through which the domain wall has not yet 
swept through. The compaction process is related to the 
growth of the dense ordered regions rather than to the 
characteristic length of the domains. The fast dynamics 
of the domain walls is thus not in contrast with the slow 
dynamics of the bulk density. 

In order to support this picture we have measured the 
fraction of persistent sites (or persistence probability) 
R{t) jl^, i.e. the fraction of sites that never changed 
their status up to time t. If the triggering process men- 
tioned above was perfectly efficient, i.e. every time a 
domain wall passes through a vacancy one triggers a pro- 
cess increasing the density locally, one would expect that 
1 — pit) is described by R{t). Otherwise the behaviour of 
l—p{t) is slower than R{t). In analogy with recent studies 
[|l4| in standard coarsening models, a very slow algebraic 
decrease of such quantities is observed as a function of 
time; we obtain in particular that R{t) (see Fig.||) scales 
as t-'^ with e ~ 0.15. The behaviour of (1 - p{t)) vs. t 
is shown in the inset of Fig.^. We observe that it is con- 
sistent with an algebraic decay 1 — p{t) ~ with an 
exponent smaller than the persistence exponent 6 which 
is compatible with our discussion above. 

Additional information on the system can be obtained 
by monitoring the following quantities. (1) A{t): activ- 
ity in the system measured as the cumulative number 
of successful moves; (2) M [t) is intended to measure the 
mobility of the domain walls. It is obtained by the cumu- 
lative sum of the absolute value of the derivative of the 
staggered magnetization in thin vertical stripes, averaged 
over all the stripes. From Figj^, which reports the results 



for both quantities, one deduces two main things: the ac- 
tivity of the system is concentrated on the domain walls 
(since A{t) and M{t) have the same functional form up to 
a constant) and both A{t) and M{t) scale as [1 — p{t))~^ . 
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FIG. 4 
for X — 0.5, 
obtained with x 
vs. t [15]. 



A(t) and AI{t) (see text for the definition) vs. t 



- 800, Ly — 50. Similar results have been 
: 0.1 and X = 0.9. The inset reports (1— p(f)) 



In order to get a better insight into the above men- 
tioned phenomenology it would be very interesting to 
have a quantitative understanding of the link between 
the coarsening process and the very slow global den- 
sity relaxation. Under the hypothesis that the system 
is translationally invariant in the direction of gravity, the 
medium is described as a one dimensional (y-averaged) 
density profile, p[x,t)^ in which particles (the domain 
wall) move. To comply with our previous description, we 
assume that the density p is only susceptible to increase 
at the positions occupied by the walkers, and remains 
quenched elsewhere. Moreover, in order to follow more 
closely what happens in the Tetris model we consider a 
situation where the motion of the walkers is coupled to 
the environment, i.e. the local density via a potential 
field depending only on the density. The problem can 
be cast in terms of two coupled equations: one describ- 
ing the over-damped motion of the walkers and another 
describing the evolution of the local density as: 



dX 
~dt 

dpix,t) 
dt 



dp{x,t) 



dp{x,t) 



--X{t) 



dx 



--X{t) 



= f{p{x,t)) 5{x~X{t)) 



Tit) 
(2) 



where r(t) is an uncorrelated Langevin force with 
{Tit)) = and (r(t)r(i')) = q6it - t'). The potential 
V attracts the walker to regions where activity has been 
intense, and repels it from unvisited regions. 

A detailed treatment of the above defined equations is 
presented elsewhere p^ . Here we summarize the main 
features. The function / should be such that p ap- 
proaches unity for long times, and hence /(p) — > as 
p ^ 1. An example of such an / is /(p) = (1 — p)", where 



3 



a is a positive exponent. The potential V should provide 
a drift toward high density regions. A simple suitable 
form is V{p) — — p"*". In fact it is possible to show that 
the latter functional form is inessential, provided V{p) 
behaves linearly in p as p — > 1. In this way ^gp^^^f*'*^ 
becomes a constant, i.e. unimportant, as p tends to 1. 
With these definitions eqs.(H) can be recast in the form 

'^^ -1-^-- ■ m, ^ = 5{x~X{€)) (3) 



F 



VF- 



dt ' ' dt 

where b ~ l/(a — 1) and we have introduced the function 
F = [l/(a— !)](! — p)^^" which, according to Eq. ^, rep- 
resents the cumulative activity on the site x (in general 
dF/dp = f{p))- In this way the choice of the function / 
is consistent with the results obtained for A{t) and M (t) 
in the Tetris model. The main idea behind this kind of 
modeling is that the high density regions (i.e. the poten- 
tial wells) tend to trap the walkers, that, in their turn 
are able to change the environment, i.e. the local den- 
sity, though their efficiency decreases with the increase 
of the density. From the combination of these two effects 
a drastic slowing down is expected. The way the walkers 
escape from the potential wells is to progressively carve 
their way out by pushing the potential barrier and so 
enlarging the compactified region. 

Different aspects come into play in the compaction pro- 
cess. One of them is related to the fact that initially 
a large number of walkers is to be introduced. Since 
our modeling only concerns the regime where a one di- 
mensional description is adapted, the typical distance be- 
tween walkers is proportional to Ly. However, when two 
walkers meet, they annihilate, and thus the subsequent 
increase the density becomes less and less effective. The 
quantification of this effect has been done in a previous 
section, through the pair correlation function of the do- 
mains. By itself, this single aspect is not sufficient to 
account for the slow densification observed numerically. 
The second aspect concerns the densification due to a sin- 
gle walker. Starting from a low density for the medium, 
we observe that the density does not remains uniform. 
Starting from any site, at low temperature, the domain 
walls first drill a potential well where they lie. How- 
ever, as the density approaches unity, the densification 
becomes less and less efficient. The only option for the 
walkers is to expand in lateral size through a progressive 
translation of the well boundaries. Specific solutions of 
this regime can be obtained as solitary waves [p^ . In 
this second regime, the densification rate is controlled by 
the velocity of the latters. Finally, the wells tend to co- 
alesce and the mean density decays as 1 — p ~ 
Though the variety of the different phenomena involved 
in the density evolution, (number of walkers, width and 
depth of potential wells, late stages crossover phenom- 
ena) renders difficult the identification of the F function 
in Eq.3, we retain that the main features observed in the 
Tetris model are captured by this walkers modeling. 



It is also important to stress that the equivalence ob- 
served between the activity and the motion of the domain 
boundaries, implies that the same treatment could also 
be carried out for particles with random shapes |^ . Here 
the existence of domains is no longer evident but the 
activity still remains confined and it is not spread out 
uniformly over the system. Finally it is worth to remark 
as our analysis could be easily exported in experimental 
setup like the one proposed in Q. Acknowledgments 
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